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] We explore, using a Ginzburg-Landau expansion of the free energy, the Larkin-Ovchinnikov-Fuldc- 

, Ferrell (LOFF) phase of QCD with three flavors, using the NJL four-fermion coupling to mimic gluon 

■ interactions. We find that, below the point where the QCD homogeneous superconductive phases 

' should give way to the normal phase, Cooper condensation of the pairs u — s and d — u is possible, 

. , but in the form of the inhomogeneous LOFF pairing. 



X 



PACS numbers: 12.38.Aw, 12.38.Lg 



I. INTRODUCTION 



| At high quark density and small temperatures Quantum-Chromo-Dynamics (QCD) predicts Cooper pairing of 
r ■ quarks due to the existence of an attractive quark interaction in the color antisymmetric channel, see Q,U Q and 
for reviews 0, E|. At extreme densities the energetically favored phase is the Color-Flavor-Locking (CFL) phase, 
characterized by a spin diquark condensate antisymmetric in both color and flavor |(| ; at intermediate densities the 
situation is much more involved, because one cannot neglect the strange quark mass and the differences 5[i in the 
ly-j , quark chemical potentials induced by (3 equilibrium. Several ground states have been considered in the literature, 
• from the 2SC phase Q , to the gaplcss phases g2SC_Q and gCFL 0, El. The gapless phases are instable as shown 
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by imaginary gluon Meissner masses (for g2SC see [lOj, for gCFL see [ill and |12|). This seems to be connected to 
the existence of gaplcss modes in these phases ^3|- An instability is present also in the 2SC phase Though 
this phase has no gapless mode, imaginary gluon masses are present when the gap A and S/j, satisfy the condition 
Qh! A/v^2 < 5 (j, < A. 

Another superconductive state discussed in the literature is the Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) Q 
phase. The relevance of this phase is based on the possibility that, for appropriate values of Sji, it can be advantageous 
for quarks to form pairs with non- vanishing total momentum: pi + P2 = 2q ^ 0, see 0] and for a review [l7j . 
As far as instability is concerned, the authors in 0] have shown that, with two flavors, the instability of 2SC implies 
that the LOFF phase is energetically favored. Moreover, in the LOFF phase with two flavors the gluon Meissner 
' masses are real [l9f . 

Thus far only the case of two species for the LOFF phase has been studied. This is not justified in QCD. At 
intermediate densities all the three quarks: u, d and s should be considered. The three flavor problem is however 
much more involved and difficult to work out. We present here a first attempt to study the three flavor LOFF phase 
of QCD. Our approach is based on a Ginzburg-Landau (GL) expansion of the free energy. Differently from the CFL 
phase, where quark matter is in (3 equilibrium while being also electrically and color neutral, here we should impose 
these conditions. We shall consider in the sequel only /3-equilibrated and electrically neutral quark matter, while 
assuming that the color-chemical potentials vanish. This is an approximation we discuss below. 



II. GAP EQUATION 

To get the gap equation in the Ginzburg Landau approximation, we start with the Lagrangean density for three 
flavor ungapped quarks: 



C 



^(ipf-M^ + ^-ro)^. (l) 

M?P = S al3 diag(0, 0, M s ) is the mass matrix and Df? = dS af3 Sij + igA a T"^Sij; /U^g is a diagonal color-flavor matrix 
depending in general on /i (the average quark chemical potential), /i e (the electron chemical potential), and fig, 
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related to color [8|. We do not require color neutrality and we work in the approximation Ha = Us = 0, which might 
be justified by the results of Ref. Q for the gCFL phase showing that yu, 3 and /Kg assume in general small values (at 
least in the region of interest, see later). Therefore in this paper 

/4f = Wii - MeQ«)^ = Mi M Q/3 > (2) 

where Q is the quark electric-charge matrix. 

We treat the strange quark mass at the leading order in the 1/fi expansion; this corresponds to a shift in the 

M 2 

chemical potential of the s quark: /i s — > fj, a — . This is the same approximation used in Refs. for the study 

of the gCFL phase. Therefore: 

2 1 1 M s 

M« = M - gMe , = A* + ^Me , Ms = A* + gMe ~ • ( 3 ) 

Another approximation we employ is the High Density Effective Theory (HDET), see [23, EI H3] and, for a 
review, Here one decomposes the quark momentum into a large component, proportional to /x, and a residual 
small component: p = /m + £; n is a unit vector and I is the small residual momentum. Moreover one introduces 
n-dependent fields ipn and ^ n by the Fourier decomposition 

^{x)= j^^{^{x) + ^ n {x)) ; (4) 

rp n and \& n correspond to positive and negative energy solutions of the Dirac equation. 
Substituting the expression (@J in the Eq. (JJJ one gets at the leading order in 1/fj, 
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^(iV-D^ + fi^Sij)^, (5) 



where = (1, n), V** = (1, — n) and pi = \ii — fi. 

It is convenient to change the basis for the spinor fields by defining ip A = {''Pun ipdg, ipbsiipdr, ipug,tpsr) il>uh, ipsn,il>db)- 
This change of basis is performed by unitary matrices Fa , whose explicit expression can be found in Ref. [ll| • To the 
Lagrangean in Eq. JSJ we add a Nambu-Jona Lasinio four fermion coupling treated in the mean field approximation. 
This corresponds to the same coupling and the same approximation used in Ref. [3]. The gap term in the resulting 
Lagrangean is conveniently treated by introducing the Nambu-Gorkov field 



so that the Lagrangean reads 



[dn fdEdt t f (E-^ + fj A )5AB -A AB (r) \ 
2^J 4tt J (2tt)2 Xa { -A* AB (r) (E + £ — fl A ) Sab J [<) 

where E is the energy, £ = I ■ n is the component of the residual momentum along n and satisfies: |£| < S, with S an 
ultraviolet cutoff. Moreover (fi) A = (fl u , p, d , p, s , p, d , p, u , p, s , Jx u , JL S , fid)- 
We assume the pairing ansatz 

3 

< i>ia c lb v« >= eQ/3/ e ^ ( 8 ) 
i=\ 

with 

A/(r) =A / exp(2iq I -r) . (9) 

In other words, for each inhomogeneous pairing we assume a Fulde-Ferrell ansatz; 2qi represents the momentum of 
the Cooper pair. The gap matrix Aab in can be expressed in terms of the three independent functions Ai(r), 
A^fr), A,s(r) describing respectively d — s, u — s and u — d pairing. The explicit expression of Aab can be found in 

MM- 
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To write down the gap equation it is useful to introduce the following components of the free quark propagator 

The quark propagator is the matrix 

whose components satisfy the Gorkov equations 

Sn = S U + S n A(r)S 21 , S 21 = S 22 A*(r)S n . (12) 

S21 is the anomalous propagator involved in the gap equation. 

The wave vectors qi should be derived by minimizing the free energy. We will fix the norms |qj| by a minimization 
procedure. As to their directions, we will limit the analysis to four structures, choosing among them the one with 
the smallest value of the energy. The first structure has all qi along the positive z— axis. The structures 2, 3, 4 have, 
respectively, qi, q2, q3 along the positive z— axis (the remaining two momenta along the negative z~ axis). This is 
obviously a limitation. It is justified by our final results that show the existence of a range of values of the strange 
quark mass where the LOFF phase, even with these limitations, is favored in comparison with other QCD phases. 

The gap equation in the HDET formalism can be written as follows Q 

A^ B (r) = i 3G V^V" ]T h* AaC h DbB j ^ j ^± j g S 21 (E,£) CD g^5 ab , (13) 

where £21 is given in Eq. I|12(l : in the above equation hr>bB is a Clcbsch-Gordan coefficient. It is expressed by the 
formula hobB = ^[F^T^Fb] in terms of the unitary matrices Fa used to write the quark fields as in ©, i.e. in the 
basis A = 1, • ■ • , 9. G is the Nambu-Jona Lasinio coupling constant, of dimension mass -2 . In what follows, we shall 
get rid of G introducing the value of the CFL gap parameter Ao as a measure of the strength of the interaction (see 
below, Eq. JTBJl). 

III. GINZBURG-LANDAU EXPANSION 

Performing the Ginzburg-Landau expansion of the propagator 

S 2 i = S 22 A*Sq 1 + Sq 2 A*S1 1 AS$ 2 A*Sq 1 + 0(A 5 ) (14) 

we get 

A/ = 11/ A/ + ^ J;j A; A 2 ; + 0(A 5 ) , 7=1,2,3. (15) 
./ 

Let us comment on the functions 11/ and Jjj appearing in this expansion. 11/ are defined as follows: IT = 
U(qi,S^ ds ) , il 2 = U(q 2 ,SiJ, us ) , II3 = U(q 3 ,S^ ud ) , with 

, fid-pu Me j- _fLs-pu Me M 2 fi e - fl d M 2 

S^ d =—^— = T , (S/itis = ^_ = ___, 5fXds = __ = __. (16) 



and 

H(q,6fj,) = 1 



2Gfi 2 





<l 


+ 6n 


log 


4{q 2 - 6n 2 ) 


) 




'1 


— <5/i 


2 8 


A 2 
^0 





(17) 



We note that II is analogous to the function determining the behavior of the free energy in the GL approximation of 
the LOFF phase with two flavors. We have introduced the parameter Ao to get rid of the ultraviolet cutoff <5. It is 
defined by 



A °^ 2SeXp ^-2G^? ■ 
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Ao is equal to the CFL gap for M s = and /i e = in the weak coupling limit, with no sextet condensation. As for 
J/j, we have, for the diagonal components: J n = J x = J{qi,5^ds), J22 = h = J(?2,<fyO> J33 = J3 = J{qz,5^ ud ) } 
with 

The off-diagonal term J12 is 

J 12 = ^!/^ I ; (20 ) 

7T 2 J 4?r (2qi -n + iJ, a - (id- ie) (2q 2 ■ n + fi s - ^ u - ie) 
J13 is obtained from J12 in (|2U[) by the exchange q 2 — > q3 and /Lt s *-> /x<j; J23 from J12 by qi — > q3 and fj, s *-* [i u . 

IV. FREE ENERGY 

Let us now consider the free energy Q. It is obtained by integrating the gap equation. The result is 
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fi = f2„ + ^(f A? + |a* + £^A 2 A 2 7 ] + 0(A 6 ) (21) 

j=i V j^/ 



with 



°- = -1^04 + ^ + ^) -1^5 (22) 

where the chemical potentials for quarks are defined in Eq. J2J and the coefficients are given by 

2(1-11,) 2 Jj 2J U 

= £ ' P 1 = Q~ ' P" = Q~ ■ ( 23 ) 

Electric neutrality is obtained by imposing the condition 



(24) 



which, together with the gap equations, gives, for each value of the strange quark mass, the electron chemical potential 
He and the gap parameters A/. Moreover one should determine qi by searching for the energetically favored solution. 
This is a complex task as it would require the simultaneous solution of the previous equations (|24|l and (|15fl together 
with: 



dQ A,^+A,5>^, 7=1,2,3. (25) 
oqi j-^ aqi 



dqi dq! f-| dq 



Moreover one should look for the most energetically favored orientations of the three vectors qi in space. A complete 
analysis is postponed to a future paper; as discussed above we have limited the analysis to the four structures 
characterized by all vectors qi parallel or antiparallcl to the same axis. Even with this limitation we are able to prove 
that there exists a window of values of M s where the LOFF phase is favored in comparison with other phases of QCD, 
as it will be seen below. As to the norms |qi|, since we work in the GL approximation, we can neglect the C(A 2 ) 

da 

terms in (12511 . As a consequence we simply get — — = 0, which, being identical to the condition for two flavors, gives 

1—1 1—1. "Ii 
the result qi = 1.1997|fytj| M,uM 

V. RESULTS AND DISCUSSION 

Our results are summarized in Figs. Ill II In Fig. Q]we give SIloff — ^norm (in units 10 6 MeV 4 ) as a function of 
M 2 //z (in MeV). 
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We report only the energetically favored solution. It corresponds to Ai = 0, A2 = A3 and q2, q3 parallel (structure 
1, which is identical in this case to the structure 2). We do not report the other solutions corresponding to local 
minima of f2. We also found a solution with Ai 7^ and A2 7^ A3 but it can only exist in a kincmatical range where 
the LOFF phase is energetically disfavored in comparison with the CFL or the gCFL phases. The structures 3 and 4 
(Q2, c[3 antiparallel) have almost vanishing gaps and correspond only to local minima of the free energy. The results 
in this figure and in the subsequent ones are obtained for fi = 500 MeV (for /j, = 400 MeV the results are qualitatively 
similar). The value of the CFL gap for zero strange quark mass is fixed at Ao = 25 MeV. This is the same value 
used in Ref. This choice, as well as the same form of the NJL coupling, with the same approximation, allows a 
comparison between our results and those of Ref. 0, see the discussion below. 




20 40 60 80 100 120 140 



M s /|Ll[MeV] 



FIG. 1: Free energy difference Qloff — fonorm in units of 10 6 MeV 4 plotted versus M^/fA (in MeV). The result is obtained 
for n = 500 MeV and Ao = 25 MeV. The line corresponds to the structures 1 and 2 with Ai = 0, A2 = A3. 

In Fig. |21 we give the gaps A7/A0 as functions of Mf//i (in MeV). The line represents the solution A2 = A3 for 
the structures 1 and 2 (in this case Ai =0). 




40 60 80 100 120 140 



M s /(I[MeV] 



FIG. 2: Gaps Aj/Ao as functions of M s 2 / fi (in MeV) for the structures 1 and 2; the curve represents Aj = A2 = A3 whereas 
Ai = 0. 



In Fig. |3 we present results for the electron chemical potential [i e . The line correspond to the energetically favored 
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solution (structure 1=2; we have neglected terms suppressed in the expansion, consistently with the HDET 
scheme) . A glance at eq. (fR))l shows that pi e is given by u e « M 2 / (4/i) , which corresponds to a symmetric splitting 
of the s and d Fermi surfaces around the u Fermi surface. Therefore in this kinematical region we have us and du 
pairings, with 

Pu + Ps = 2 q 2 , Pu + Pd = 2 q 3 = 2 q 2 . 

The gaps A2 {us pairing) and A3 {ud pairing) have to be almost equal since they depend only on the absolute values 
of the splittings, see Fig. Since the separation between the d and s Fermi surfaces is higher, one does not expect 
ds pairing, which is confirmed by the result Ai = in the region where LOFF prevails. 

H e [MeV] 

40 
30 
20 
10 

20 40 60 80 100120140 

M s /|i [MeV] 



FIG. 3: The electron chemical potential fi e as a function of M^/u. Units are MeV. The line corresponds to A2 = A3, Ai = 0. 



In Fig. 0]wc present comparison of different phases of QCD. In order to comment this figure, let us start assuming that 
all the other phases are stable, meaning that in some way it is possible to cure the instability due to the imaginary 
gluon masses. In this case, following the graph for decreasing values of M s 2 //j, we see that at about M^/n = 150 
MeV the LOFF phase has a free energy lower than the normal one. This is a second order transition as it can be 
seen from Fig. Then the LOFF state is energetically favored till the point where it meets the gCFL line at about 
Mil a = 128 MeV. This is a first order transition since all the gaps are different in the two phases (for the gCFL case, 
see [8j). Then the system stays in the gCFL phase up to M 2 //x « 48 MeV where it turns into the CFL phase via a 
second order transition (see Q)- 

However, if the gapless phases are unstable, then they should not be considered, and the LOFF phase is the stable 
phase from M%/(j, = 150 MeV down to about M 2 //i = 90 MeV where the LOFF line meets the CFL line, with a first 
order transition (this can be seen by comparing our gaps with the Ac fl ~ 23 MeV at this value of M 2 / fj,) . 

We should also add that at the moment it is still unknown if the LOFF phase with three flavors suffers of chromo- 
magnetic instabilities. This problem is left to future investigations. 

VI. CONCLUSION AND OUTLOOK 

We have explored in the framework of the Ginzburg-Landau expansion the LOFF phase of QCD with three flavors, 
using the NJL four-fermion coupling to mimic the gluon interactions. We have worked on the ansatz of a single 
plane wave behaviour for each quark pairing, which is the simplest generalization of the gCFL phase that takes into 
account the possibility of anisotropic condensation. We found that near the point where the CFL phase should give 
way to the normal phase, Cooper condensation takes place in the form of the LOFF pairing. Our analysis has some 
limitations. First, we have assumed vanishing color chemical potentials /13, /is- The results of Ref. |9J show that 
in the region where the LOFF state dominates the color chemical potential have rather small values, in particular 
smaller than /i e . However non vanishing values of fj.3 and /is are expected to increase the LOFF free energy and 
therefore a more complete calculation is necessary. Second, we have considered the three possible momenta qi along 
the same direction. Third, more than one plane wave might be present in each condensate. Finally we have treated 
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FIG. 4: Free energy differences Qloff — Qnorm in units of 10 6 MeV 4 plotted versus Ml / fj, (in MeV) for various QCD phases. 

the strange quark mass at its leading effect, i.e. by a shift in its chemical potential, which is also an approximation 
[22I ]. We plan to address all these issues by a more refined study in the future. 
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